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We present an experimental study of an optomeclianical system formed by a vibrating tiiin semi- 
transparent membrane within a high-finesse optical cavity. We show that the coupling between 
the optical cavity modes and the vibrational modes of the membrane can be tuned by varying 
the membrane position and orientation. In particular we demonstrate a large quadratic dispersive 
optomechanical coupling in correspondence with avoided crossings between optical cavity modes 
weakly coupled by scattering at the membrane surface. The experimental results are well explained 
by a first order perturbation treatment of the cavity eigenmodes. 
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I. INTRODUCTION 

The study of cavity optomechanics has recently 
sparkled the interest of a broad scientific community 
due to its different applications, ranging from sensing of 
masses, forces and displacements at the ultimate quan- 
tum limits [U [2] , to the realization of quantum interfaces 
for quantum information networks 3-6 , up to tests of 
the validity of quantum mechanics at macroscopic level 
[3 IB] • A large variety of devices has been recently pro- 
posed and tested in which a driven cavity mode interacts 
with a mechanical resonator due to the fact that the cav- 
ity mode frequency w depends upon the effective posi- 
tion of the mechanical element z [31 [51-[T5] . In most cases 
one has a linear dependence between w and z, so that 
the optical field exerts an homogeneous force on the res- 
onator, associated with either radiation pressure or the 
gradient dipole force, which has been proposed and used 
for cooling the resonator motion [I4l420j . In this case 
the phase shift of the output light is proportional to the 
mechanical displacement and one can implement high- 
sensitive readout of forces and displacements [5T1 [52] . In 
the so-called membrane-in-the-middle scheme of Ref. [12] 
the position dependence of the cavity mode frequency 
is caused by a semi-transparent thin membrane placed 
within the cavity. In such a case, when the membrane 
is placed at a node or at an antinode of the cavity field, 
uj{z) is quadratic in the position of the mechanical ele- 
ment, and one has a dispersive interaction, which allows 
for new nonlinear optomechanical phenomena [23], and 
for the quantum non-demolition detection of the vibra- 
tional quanta of the membrane motion [T^J [25] . 

Here we present an experimental study of the 
membrane-in-the-middle scheme. We show how the op- 
tomechanical coupling can be tuned by varying the po- 
sition and the orientation of the membrane with re- 
spect to the cavity axis. In particular, we show that 
by appropriately tilting the membrane one can couple 
the various Hermite-Gauss modes of the cavity and in- 
duce avoided crossings between the frequency shifts of 



the cavity modes (see also Ref. [25 ). When the split- 
ting at an avoided crossing is very small, the correspond- 
ing value of the second-order derivative of the frequency 
with respect to the position d'^ijj{z) / dz"^ becomes very 
large and one achieves a significative quadratic disper- 
sive optomechanical interaction. Here we experimentally 
achieve d^uilz) / dz^ — 27rx4.46 MHz/nm^, which is com- 
parable to the value achieved in Ref. [5S]. The paper 
will also explain the experimental results by means of a 
first order perturbation theory able to illustrate how the 
tilted membrane determines the new cavity eigenmodes 
and their frequency shifts. 

The outline of the paper is as follows. In Sec. II we il- 
lustrate the main aspects of the membrane-in-the-middle 
scheme, while in Sec. HI we present our experimental 
setup. In Sec. IV we present the experimental results 
on the tunability of the optomechanical coupling and the 
realization of a strong quadratic coupling, while Sec. V 
is for concluding remarks. In the Appendix we provide 
details on the perturbation theory able to explain quali- 
tatively and quantitatively the modification of the cavity 
modes caused by the membrane. 

II. THE MEMBRANE-IN-THE-MIDDLE 
SCHEME 

The membrane-in-the-middle system is formed by a 
Fabry-Perot optical cavity with a thin semi-transparent 
membrane inside. The empty cavity supports an infi- 
nite set of optical modes, conveniently described by the 
Hermite-Gauss modes [26]. The thin membrane is a di- 
electric slab of thickness Lj^ and complex index of re- 
fraction rijv/ = rifl -|- in/. When it is placed within the 
cavity at a position zq measured along the cavity axis, 
the mode functions and their frequency change in a way 
which is dependent upon the position and orientation of 
the membrane with respect to the cavity [T21 [5S1 [571 HH] ■ 

Considering the membrane motion means assuming 
that its mean position along the cavity axis oscillates 
in time, zq — + z{x,y,t), where z(x,y,t) gives the 
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membrane transverse deformation field, and is given by a 
superposition of the vibrational normal modes. One typ- 
ically assumes the high stress regime of a taut membrane, 
in which bending effects are negligible and the classical 
wave equation well describes these normal modes 28.,.29j. 
The membrane vibrational modes are coupled to the opti- 
cal cavity modes by radiation pressure, and therefore one 
has in general a multimode bosonic system in which many 
mechanical and optical modes interact in a nonlinear way. 
However, one can often adopt a simplified description 
based on a single cavity mode interacting with a single 
mechanical mode [HI UHl 130] • This is possible when: i) 
the driving laser mainly populates a single cavity mode 
(with annihilation operator a), and scattering into other 
modes is negligible [21]); h) the detection bandwidth is 
chosen so that it includes only a single, isolated, mechan- 
ical resonance with frequency flm (described by dimen- 
sionless position q and momentum p operators, such that 
= i). By explicitly including cavity driving by a 
laser with frequency lul and input power V, one ends up 
with the following cavity optomechanical Hamiltonian, 
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where E — ^IVKi^jhjjL^ with kq the cavity mode band- 
width in the absence of the membrane. In Eq. ([T]) we 
have included the radiation pressure interaction within 
the cavity mode energy term, by introducing a position- 
dependent cavity frequency u;(g), which can be written 
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where loq is the cavity mode frequency in the absence 
of the membrane, and [2:0(9), a^,, ay] is the frequency 
shift caused by the insertion of the membrane. This 
shift depends upon the membrane position along the 
cavity axis 2:o('?), which in turn depends upon the co- 
ordinate q because 2:0(9) = 20 -f xoOg, where zo is the 
membrane center-of-mass position along the cavity axis, 
= \/hJrn£L^n is the spatial width of the mechanical 
zero point motion (m is the effective mass of the mechan- 
ical mode), and is the dimensionless overlap integral 
between the transverse mode functions of the selected 
mechanical and optical modes [55]. The frequency shift 
also depends upon ax and cty, the tilting angles around 
the X and y axis respectively. The (x^y^z) axes form a 
left-handed cartesian frame with the origin at the center 
of the cavity, and x and y are otherwise arbitrary, due to 
the cylindrical symmetry of the optomechanical system 
around the cavity axis z. 

When the thin membrane is perfectly aligned, = 
ay = 0, and is placed very close to the cavity waist, 
the Hermite-Gauss modes still represent the cavity eigen- 
modes with a very good approximation, because their 
wavefronts fit well with the membrane. The frequency 
shift in this case is a simple periodic function of 2:0, 
which is maximum at the antinodes and minimum at 



the nodes of the intracavity field, and mode degeneracy 
is not removed [28l |30]. When the membrane is ap- 
preciably shifted from the waist and/or is tilted, light 
scattering of the Hermite-Gauss modes at the membrane 
surface is no more negligible and the cavity eigenmodes 
are significantly modified [25]. When the longitudinal 
shift from the waist and the tilting angles are sufficiently 
small, one can describe the situation by means of a degen- 
erate first-order perturbation theory of the wave equa- 
tion within the cavity. In this perturbation limit, the 
new eigenmodes are linear combinations of few Hermite- 
Gauss modes and the corresponding frequency shifts 
Suj [zo{q),ax,ay] can be correspondingly evaluated [32j . 
In this paper we will study both experimentally and the- 
oretically the behavior of these frequency shifts: we will 
see in particular that the cavity mode frequencies, and 
therefore the optomechanical coupling as well, can be 
fine-tuned and controlled in situ by varying the position 
and the orientation of the membrane. 

In fact, the optomechanical coupling between the se- 
lected cavity and membrane vibrational modes is pro- 
vided by the first-order (and eventually higher-order) 
term in the expansion of Suj [20(9): ctx, cty] as a function 
of q. At most membrane positions zq a first order ex- 
pansion in q of Eq. ^ provides an accurate description 
of the physics: in this case one has a standard radiation 
pressure optomechanical interaction with a single-photon 
optomechanical coupling strength given by 
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where duj/dzo can be directly measured experimentally, 
while xq and Q depend upon the chosen optical and me- 
chanical modes. 

Instead, when the membrane center 20 is placed ex- 
actly at a node or at an antinode of the cavity field, 
or at an avoided crossing between nearby frequencies, 
the first-order term in the expansion of uj{q) vanishes, 
and one has to consider the higher-order term, which 
is quadratic in q. This latter term describes a disper- 
sive interaction between the optical and the vibrational 
modes whose coupling rate is given by the second or- 
der derivative d'^uj{q) / dcf' . This unique property of the 
membrane-in-the-middle scheme has been first discussed 
in [T2I [3^ and if such a quadratic coupling is sufhciently 
strong, could be exploited for a quantum non-demolition 
measurement of the vibrational energy |121 124j . 



III. THE EXPERIMENTAL SETUP 

Our membrane-in-the-middle setup is schematically 
described in Fig. I. Laser light is produced by a Nd:YAG 
laser (Innolight) with wavelength A — 1064 nm. The light 
passes through an optical isolator, and is sent to the op- 
tomechanical cavity via two steering mirrors. The cavity 
is formed by two dielectric mirrors each with a radius 
of curvature Ri = R2 = 10 cm (coated by Advanced 
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Thin Films) and separated by a distance L = 9 cm. The 
measured cavity finesse without the membrane is equal to 
Fq ~ 66000. The membrane is mounted on a piezo-motor 
driven optical mount that controls the angular alignment 
of the membrane. The optical mount is in turn mounted 
on xyz stack of piezo-motor driven linear stages that are 
used for moving the membrane in space. The stages are 
vacuum compatible, two of them are used for centering 
the membrane with respect to the optical axis, while the 
remaining one positions the membrane along it. Rota- 
tion around two axes perpendicular to the optical one 
provides us with tip and tilt control. 
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FIG. 1: Schematic description of the experimental setup. 

The membrane used in the experiment is a commer- 
cial, 1 mm X 1 mm Si3N4 stoichiometric x-ray window 
(Norcada), with nominal thickness = 50 nm, and in- 
dex of refraction nn ~ 2. The membrane is supported 
by a 200 fim thick Si frame, and it has been chosen due 
to its high mechanical quality factor and very low optical 
absorption (n/ ~ 10~^ at A = 1064 nm), as discussed 
in Ref. |33) . This parameter corresponds to an inten- 
sity reflection coefficient TZ — 0.18, which has been also 
experimentally verified. 

The Fabry-Perot optical cavity is mounted inside a 
custom made vacuum chamber which hosts also the mem- 
brane. The chamber is pumped down with a turbo- 
molecular pump down to 10^^ mbar. Once the base 
pressure is reached, the pumping is switched off and 
the chamber is disconnected from the pumping station, 
and the system is left in static vacuum. Usually during 
the measurements the static vacuum inside the chamber 
reaches values as high as 10^^ mbar. 



IV. MEASUREMENT OF THE CAVITY MODE 
FREQUENCY SHIFTS 

We measure the optical cavity modes excited by the 
driving laser by monitoring the optical power transmit- 
ted through the cavity using the photodiode shown in 
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FIG. 2: Predicted frequency shift Ai/ — 5w/27r of the 
Hermite-Gauss modes close to the selected TEMoo,p (blue 
curve, p = longitudinal number) vs the membrane posi- 
tion along the cavity axis, for a perfectly aligned membrane 
around the waist. The green curve refers to the TEM2o,p-i, 
TEMii,p_i, TEMii,p_i degenerate triplet, and the violet 
curve to the TEM4o,p-2, TEM31.P-2, TEM22,p-2, TEMi3,p-2, 
TEMo4,p_2 degenerate quintuplet. See Sec. |III| for the other 
system parameters. 



Fig. 1. We can also determine the transverse profile of 
the cavity mode by imaging the transmitted beam with 
a video camera. The laser is tuned and mode-matched 
to one of the cavity's TEMqo mode, and its frequency 
is then scanned over a prefixed range. In such a case, 
and if the input power is not too low, also higher order 
transverse TEM,„„ modes (with not too large m and n), 
which are close in frequency to the driven TEMqo rnode, 
are populated and can be detected. For the system pa- 



rameters given in Sec. Ill the relevant cavity modes are 
shown in Fig. 2, which refers to a perfectly aligned mem- 
brane = Oiy = 0, whose position is scanned for one 
half of the laser wavelength around the cavity waist. 

As discussed above, if the membrane is slightly mis- 
aligned, ax, ay 0, the Hermite-Gauss TEMm„ modes 
are no more cavity eigenmodes, and the frequency 
shifts are consequently modified. In particular: i) 
degeneracies are removed; ii) avoided crossings appear 
in correspondence of cavity modes which are coupled by 
the perturbation caused by the tilted membrane. This 
fact has been experimentally verified and it is shown 
in Fig. 3, where the optical power transmitted through 
the cavity as a function of the laser frequency and the 
membrane position is shown. Fig. 3 refers to a 0.5 
mm shift from the cavity waist and to tilting angles 
ax — —0.21 mrad, ay = 0.15 mrad. The crossings of the 
TEM2o,p-i triplet with the TEMoo,p mode and with the 
TEM4o,p_2 quintuplet are visible, together with the re- 
moved degeneracies (p denotes the longitudinal number 
of the driven TEMqo mode) . In the zoomed part avoided 
crossings are visible and their behavior is satisfactorily 
reproduced by the red curves, which correspond to the 
theoretical prediction of the perturbation theory devel- 
oped in the appendix. In particular, the upper avoiding 
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fact that cylindrical symmetry around the cavity axis 
is satisfied at a very good level of approximation by 
our setup. We have checked this fact experimentally, 
by verifying that one gets essentially the same results 
if the values of ax and ay are exchanged. This means 
that in practice, for not too large misalignments, only 
a single tilting angle around an appropriate transverse 
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Fig. 3 shows that, as expected, for most membrane 
position along the cavity axis, the frequency shift = 
(5w/27r of the various optical modes has a linear de- 
pendence on zo and therefore on q, corresponding to 
the usual radiation pressure optomechanical interaction. 
Fig. 3 also shows that the largest coupling strength is 
achieved for the TEMoo.p mode for zg values halfway 
between a node and an antinode, and it is given by 
\dijj/dz[j\ = 27r X 2.8 MHz/nm, which essentially coin- 
cides with the maximum achievable value \du) / dzQ\max = 
2ui{^\/TZ/L (see Ref. The membrane we are cur- 

rently investigating has a vibrational mode with effective 
mass m ~ 34 ng and ~ 380 kHz, corresponding to 
a zero-point motion width xq ~ 1.11 x 10~^ nm which, 
assuming an optimized overlap integral ~ 1, implies a 
single-photon optomechanical coupling G'o/27r ~ 3.3 Hz. 

At membrane positions zq exactly corresponding to 
nodes, antinodes and avoided crossings, the linear term 
in q is zero and one has a purely quadratic dispersive 
optomechanical interaction. 
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FIG. 3: Optical power transmitted through the cavity as a 
function of the laser frequency and the membrane position, 
in the case of membrane tilting angle = —0.21 mrad, 
ay = 0.15 mrad, and shifted by 0.5 mm from the waist. The 
crossings of the TEM2o,p-i triplet with the TEMoo,p mode 
and with the TEM4o,p-2 quintuplet are visible, together with 
the removed degeneracies and the presence of avoided cross- 
ings. Red curves in the zoomed part refer to the theoretical 
prediction of the perturbation theory developed in the Ap- 
pendix. 



crossing concerns the TEMoo,p and TEM2o,p-i modes 
because the other two modes of the triplet, TEMii_p_i 
and TEMo2,p-i, are not coupled to the TEMoo.p by 
the perturbation caused by the membrane. In the 
lower part of the figure one has instead three avoided 
crossings, associated to three pairs of modes which 
are coupled by the membrane perturbation, i.e., the 
{TEM2o,p-i, TEM4o,p-2}, {TEMn,p_i, TEM3i,p_2}, 
and {TEMo2,p-i,TEMo4,p-2} pairs. In fact, the tilted 
membrane does not couple other pairs of triplet and 
quintuplet modes. This behavior is consistent with the 



-205 -175 From Fig. 3 one gets that |9^cij/9zq| ~ 27r x 24 kHz/nm^ 

at nodes and antinodes of the TEMoo.p mode, while one 
has a significant increase, by almost two orders of magni- 
tude, at the avoided crossing between the TEMoo,p mode 
and the triplet, where we measured \d'^u}/dzQ\ ~ 27r x 1.0 
MHz/nm^. This shows that avoided crossings allow to 
achieve significant values of the dispersive optomechani- 
cal coupling G2 :25j, and suggests that by appropriately 
adjusting the membrane shift from the cavity waist and 
its tilting angle, one could "engineer" the spectrum of 
the optical cavity modes and generate avoided crossings 
with a very large value of the second order derivative 
|9^a;/9zo|. In fact, |9^w/i9zo| is inversely proportional 
to the frequency splitting at the avoided crossing point, 
and therefore one should look for avoided crossings with 
a very small splitting, that is, between two cavity modes 
which are very weakly coupled by light scattering at the 
membrane surface. 

We have first experimentally investigated how the fre- 
quency splitting and the associated second order deriva- 
tive at the avoided crossing varies as a function of the tilt- 
ing angle ax, while keeping all the other parameters fixed. 
The coupling between the TEMoo.p and the TEM2o,p-i 
modes decreases for decreasing \ax\, and therefore one ex- 
pects to find a minimum of the splitting Aac, and a maxi- 
mum of the second order derivative around ax = 0. This 
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FIG. 4: Frequency splitting between the TEMoo.p and the 
TEM2o,p-i modes at the avoided crossing, Aac (a), and sec- 
ond order derivative of the cavity mode frequency (b) vs the 
tilting angle Ox- The other parameters are the same as in 
Fig. 3. In both figures the full curves refer to the theoret- 
ical prediction of the perturbation theory developed in the 
Appendix. 



behavior is confirmed by the experimental data shown 
in Fig. 4, which refer to the avoided crossing between 
the TEMoo,p and the TEM2o,p-i modes, Aac, and to the 
same parameters of Fig. 3, except that we have varied 
around Ux = 0. In particular we see that around a^; ~ 
we achieve \d'^u}/dzQ\ ~ 27r x 2.0 MHz/nm^. 

We have then looked for larger values of the quadratic 
optomechanical coupling in different parameter regions. 
We have found an interesting configuration by combin- 
ing the effect of increasing the shift of the membrane 
position from the cavity waist at the center with that 
of membrane tilting. When the membrane is apprecia- 
bly shifted from the waist, the curved wavefronts of the 
cavity modes do not fit anymore with the flat membrane 
surface, and one has an additional perturbation which, 
differently from membrane tilting, does not break the 
cylindrical symmetry around the cavity axis. The com- 
bined effects of the two perturbations may lead to very 
small couplings and therefore to avoided crossings with 
large second order derivative. Such a situation is shown 
in Fig. 5, which corresponds to tilting angles — 0.77 
mrad, ckj, ~ 0, and to a membrane shifted by 1.2 mm 



from the waist. The dots refer to the measured values of 
the frequency shifts of the appreciably populated modes. 
The prediction of the perturbation theory developed in 
the appendix (given by the full lines in Fig. 5) well repro- 
duces the frequency shifts provided that the quintuplet 
modes are included, even if they are never significantly 
populated. 
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FIG. 5: The dots refer to the measured central peak frequency 
shifts Ai/ = 5Lj/2n of the appreciably populated cavity modes 
vs the membrane position along the cavity axis. The four dif- 
ferent colors, grey, red, blue and green, correspond to the pop- 
ulated modes and can be associated with the TEMoo,p mode 
and the triplet modes, TEM2o,p-i, TEMii,p_i, TEMo2,p-i, 
only far from crossing points (see the corresponding trans- 
verse mode images). The full lines corresponds to theoretical 
prediction of the perturbation theory developed in the ap- 
pendix for the corresponding membrane configuration, which 
is ax ~ 0.77 mrad, ~ 0, and shifted by 1.2 mm from 
the waist. The quintuplet modes are not appreciably popu- 
lated in this case, but they need to be taken into account in 
the perturbation theory in order to reproduce satisfactorily 
the frequency shifts and some of them are in fact shown in 
the plot. The avoided crossing of interest, with a very small 
splitting, is the one between the grey and red curve around 
zo = 0. 

The interesting point in Fig. 5 is the avoided crossing 
between the grey and the red curve, around zq = 0, 
which is characterized by a very small splitting. The 
zoom of this part is shown in Fig. 6 where, together with 
the measured and predicted frequency shifts, we show 
the images of the transverse pattern of the modes at var- 
ious membrane positions. The splitting at the avoided 
crossing is of the order of 1 MHz and, more important, 
one has at this point \d'^uj/dz^\ = 27r x 4.46 MHz/nm^, 
which is significantly larger than the values at nodes 
and antinodes, and is comparable to the value measured 
in Ref. ^W. In this avoided crossing point, the splitting 
is very small because the combined effect of tilting and 
shift from the waist causes a small but nonzero coupling 
between two orthogonal linear combinations of triplet 
and the TEMqo modes. These two cavity eigenmodes 
have the transverse pattern shown in Fig. 6, and with 
a good approximation they can be written as \4>)red = 
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FIG. 6: Zoom of Fig. 5 around the avoided crossing point in 
the vicinity of the position zo = 0. The measured frequency 
shifts (dots), the theoretical predictions (fuU lines) are shown, 
together with the image of the transverse pattern of the new 
cavity eigenmodes in various membrane positions. 



Appendix A: First order degenerate perturbation 
theory 

To describe the intracavity electromagnetic field in the 
presence of a slightly tilted membrane we adopt a de- 
generate first-order perturbation theory similar to that 
discussed in Ref. [32] and start from the Gauss-Hermite 
modes of the empty cavity as zero-th order solution. 
We consider a symmetric cavity formed by two identical 
spherical mirrors with radius of curvature R, separated 
by a distance L, in the coordinate system with the z axis 
along the cavity axis and centered at the cavity center. 
The thin membrane is a dielectric slab of thickness 
and complex index of refraction hm = nn + in/ . 

The time-independent wave equation of the empty in- 
tracavity electromagnetic field is 



[\TEM2o,p-i) + \TEMo2.p^i) ~ \TEM, 
" , = [2\TEMo2,p-i) - \TEM2o,p-i 



I gray 



'oo,p>] /V3, and 
F \TEMoo,p)] /V6. 



V. CONCLUSIONS 

We have studied both theoretically and experimen- 
tally a membrane-in-the-middle setup formed by a high- 
finesse Fabry-Perot cavity with a thin semi-transparent 
SiN membrane inside. We have seen that the position 
and orientation of the membrane within the cavity al- 
lows to fine-tune the frequencies of the cavity modes, and 
through it, the optomechanical coupling of these modes 
with the vibrational modes of the membrane. In most 
membrane positions the frequency shifts of the cavity 
mode depends linearly upon the membrane deformation 
and therefore one has the traditional radiation pressure 
coupling between optical and mechanical modes. How- 
ever at the nodes and antinodes of the cavity field the 
linear term vanishes and one has a dispersive interaction, 
which is quadratic in the position operator of the mechan- 
ical mode. We have demonstrated that such a quadratic 
coupling can be enhanced by two orders of magnitude 
in correspondence of avoided crossings with small fre- 
quency splitting. We have demonstrated a quadratic cou- 
pling term comparable to that achieved in [25] , associated 
with an avoided crossing between two new cavity eigen- 
modes which are linear combinations of the TEMoo,p, 
TEM2o,p-i, and TEMo2,p-i modes. 
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where we have defined 4>j{r) — \/2Yie[(f)j{r)] for future 
convenience, j is a collective index corresponding to the 
triplet of natural numbers {lj,mj,nj), which determines 
the j-th Gauss-Hermite mode [26) . 

Defining (f>Y' = ■{ * ^ ^ , we have 
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Hf„ is the m-th Hermite polynomial, 

w']p ^Jl + [z/zny gives the j-th 

transversal shape along z axis, w'^p — ^J^zjifkj is the 
j-th Gaussian beam radius at the cavity waist, with zr — 
(L/2)^(l + g)/(l — g) the Rayleigh range of the optical 
cavity, and kj = (tt/L) ^Ij + {rrij + rij + 1) arccos(5)/7r] is 
the norm of the j-th wave vector [2S] . The set of functions 

forms an approximate complete set of orthonor- 
mal functions for the space region within the Fabry-Perot 
cavity Q, also satisfying null boundary conditions at the 
surface of the end cavity mirrors, that is, 



d^r(f),{r)4)j{r) w S^j ^ Si^i.S.m.,njSninj , (A5) 
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The insertion of the membrane in a tilted and shifted 
position modifies the time-independent electromagnetic 
field wave equation in this way 



(A7) 



7 



where V{r) = (n|,j — l) rectj [2 



Zc{x,y)]/ta^^}, where 
rect(C) is defined as rect(C) = if |C| > 1/2, and 
rect(C) = 1 if Id < 1/2; moreover t^^ is the cor- 

1 /2 

rected thickness defined as ta^^ — Ld/ cos[(a^ ] • 

Zc{x, y) — zo + a,j;X + ayy gives, approximately, the z-axis 
projection of a generic point, of transversal coordinate 
{x,y), belonging to the (near orthogonal to z axis) cen- 
tral plane of symmetry of the membrane. The V function 
gives therefore the perturbation caused by the presence 
of the membrane within the cavity. 



We write the solution ip of the equation ( A7 1 as a linear 
combination of the orthonormal basis functions {<l>j} 



(A8) 



with Vij = fgd^r(j)i{r)V{r)(t>j{r), the matrix elements 
of the V function in the chosen orthonormal basis (see 
subsection |A 2| for an approximated explicit calculation 
of 

In the simpler case of two degenerate cavity modes 
with frequencies uJi and 012 and with degeneracies ni and 
respectively, the system of equations ( A9 ) is a finite 
system with n = rii + n2 equations 



,(i-fcf/fc2) + ^c,y,, 



0, i £ {1,2,. ..,n}. (AlO) 



Substituting ( A8) in ( A7), using ( Al I, multiplying by 



integrating over Q and using ( A5 1 , we obtain the equa- 
tions for each i, 



c4l-kf/k')+Y,c,V, 

j 



0, 



(A9) 



We have k'j/kf ~ 1 for j E {1, 2, . . . , ni}, and we define 
?7 = fcf /fcj = fc|/fcf for j E {ni + 1, ni -|- 2, . . . , n}, where 



typically rj fa 1. Defining A 



for j E {1,2,.. 

as 



fcy fcf an d Wjj 



1 + V, 



n}, we can write (AlO) in matrix form 



J3 



Vi, 



+ 1 



Kii+l,n 



Cni + 1 








(All) 



which is a sort of an eigenvalue problem with respect to A. Solving (All I gives A as a function of the membrane 
position specified by {zq, ax, cty). In particular we obtain the frequency shifts with respect to the frequency of the 
driving laser wq. 



5uj{zo,ax,ay) = cfci[A(zo, "a;, "y)] 



(A12) 



In the case of stoichiometric Si3N4 membrane one has a negligible optical absorption, that is nn « 2 and n/ ~ 10 ^. 
Defining the extinction coefficient k = nj/nj^, one has k < 10~^ ^ 1, and one can exploit the first order approximation 



X[nM{K)] w A[nM(0)] + {dnMHnM{K)]}n\,j{K)\^^^K = A(nfl) 4- 2miX'{nji). 



(A13) 



Using this approximation, we can calculate the A's as the 
real roots of the determinant of the matrix of coefficients 
in ( [All ) with real entries (each element of the matrix 
with um — nn). 



1. A singlet-triplet case study 

We apply now the general method exposed above to 
the particular case of a singlet mode coupled to a triplet 
of modes. That is, we assume, as it is relevant for the 
experimental system described in this paper, that the in- 



put laser drives a non-degenerate TEMoo mode with lon- 
gitudinal index I, specified by {(Z,0,0)}, and associated 
with the non degenerate frequency wi = (ttc/L) x [/-|- 
arccos (g) /tt] . This mode is coupled to the triplet speci- 
fied by the modes {(/-I, 2,0), (/-1, 1,1), (;- 1,0,2)}, 
associated with the three-fold degenerate frequency 012 = 
(ttc/L) X [(/ — 1) -I- Sarccos (5) /tt] . 

We have wi « W2 « oj^ by hypothesis. As discussed 
in the text, due to the cylindrical symmetry with re- 
spect to the cavity axis, we can consider (without loss of 
generality) one of the two (small) angles, ax and ay, to 
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be zero (we choose ay = 0). In particular, by defining 
Co = zo/zn, we consider two sub-cases: - small arbitrary 
Co s-nd a^] - small arbitrary Co with = 0. 



a. One- angle-tilted and shifted membrane 

If the membrane is generally not aligned (with a.y — 0) 
and not centered at the waist (small arbitrary Co and ax), 
by applying (All), we obtain the eigenvalue equation 
det[V„{\)] 



for A, where 

A Vi2 

W22 --qX V24 
14^33 - jyA 
V2i W44 - rjXj 

(A14) 

The solution of the eigenvalue equation gives: Ai = (l + 
^33)777 and three others values {Ai, A2, A3} which are 

given by the roots of the cubic equation Yl\=o 0,1^'^ — 0, 
where: 




Co 

ai 

0.2 



(A15) 
(A16) 



-77(771^11 + W22 

-VVl2{Vi2 + Vu 

+ W44{vWn + W22), (A17) 

Vi|(W^44 - V24) + VuVx2{W22 - V24) + 

+ Wn (F24 - W22W44) , (A18) 



Aligned and shifted membrane 



If the membrane is perfectly aligned {ax = ay = 0) 
and not centered at the cavity waist (small arbitrary Co)j 
the eigenvalue equation simplifies because only one off- 
diagonal matrix element is nonzero, corresponding to the 
(equal) couphng of the TEMqo mode with the TEM20 and 
TEM02 modes, i.e.. 



14^11 - A 


V12 





Vi2 


to 


W22 - rjX 














1^22 - IIX 





V12 








W22 - 



Vr,iX) 

\ " " I-' 22 ^ " 

7?Ay 
(A20) 

The eigenvalues can be explicitly obtained and are given 

by 



Ai 



^3,4 



A2 - 



1 + V^2 



22 



V 



± 



V 



1 + Vn 



-.1/2 



(A21) 



. (A22) 



1 + Vu, 7e{i,...,4}. 



(A19) 



. - . . . - - 

A further simplification is that in this case also the matrix 

elements Vij results less involved than the previous case 
with a nonzero tilting angle. 



2. Approximate evaluation of the matrix elements 



It is possible to derive an approximated explicit expression for the matrix elements Vij. We can write: 



(A23) 



= (nil ~ l)Re[y^dV</),(r)C(r)0,(r) + d^r 0,(r)C(r)(/.*(r)J - {nl, - l)Re[/{;-^^ (Co, a,, a,) + '^(Co, a., , 
where we have defined 



ll;\Co,ax,ay)= / d'r4'-'\r)C{r)^^f\r). 



(A24) 



The problem is now reduced to the evaluation of the integral (A24) and of its real part in particular. Using the fact 
that Co is small one gets (see also [32] ) 



is) 
ijOO 



iAK 



is) 



ij02 "t" ^ ij20)^0 -r \^ iji2 ~r i ijsojaxzj "t" ^03 ^ ij2l)^Vij ( i 



(A25) 



where, by introducing 5kij = {ki~ kj)/{ki + kj), — {ki±kj) / {2zji), we have used the notations raf^ = {mi±mj)/2, 
n% = (7i,±7.,)/2, iff;' = (l+.s)[4-(7n++n+. + l)]/2+(l-.s)[At-.-(mr.+r,r.)]/2, Aiff;' = [l + .s + (1 - s)fk^ /2, 



ax 



ax/ J K^-, ayij = ay/ J k^-, £\f = [ii -1-1- s{£j - 1)] /2, and also defined 



^ifqpi"^ ,ay) ^ J{mi,mj, K\f axij ,q,Skij)J{ni,nj, K^"^' ay,j ,p,5kij) 



(A26) 
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where 

J(m, n, c, q, d) 



+ C30 



dxe {x — ic^ Hm [Vl + S{x — ic)]H„ [Vl — S{x — ic)] = J^(m,n,c,q,S) — iJ^{m,n,c,q,6). 

(A27) 

This integral can be exphcitly evaluated using the properties of the Hermite polynomials, and one gets 

J«(m, n, c, q, 6) = <9™ar, {gs [cos(2/5c)i?e - sin{Vsc)Ro\ j^^^t.^o' (^28) 
J\m, n, c, g, <5) = a^a,", {gs [sm{2fsc)R., + cos{2fsc)R,] j.^^.^^Q, (A29) 

where fs{ti,t2) + (5 + ta^l - 35(^1,^2) = expj-i^ -tl + /|}, and 



L«/2J 



fe=0 ^ 

L(9-1)/2J 
Ro{q,fs,c) = ("1) 



fc=0 



2k 



(A30) 
(A31) 



with Skifs) ~ "121=0^ {2i)^(^ + '^/'^)fs ■ Using these results, one finally gets the following expression for the 



rk-2] 



relevant integral of Eq. (A25), 



Re[4*^(Co, ax, cty)] = xi cosx2|Re 



r. 



^ ij21 ^ ^ ij03 

p(s) , p(s) 
ijl2 ij30 



} 



xi sina;2 



< Im 



(^) 



Im 



p(s) I p(s) 
ij02 ' ij20 



Re 



ii21 "I" ^ ij03 



ijOQ 

O-yij 



AK, 



Re 



Co + Im 



pC^') 

J- ijl2 



I J 30 



ij02 ' ^ ij20 



Co 



+Im 
+Re 

where we have defined 

x^ - {7rH< +»;■ \m, !n, !) ' (t„^„ /i) sine (i^ff t^., /zj,) e-<°" - 
X2 = 2KlfCo + Tri\f. 



(A32) 

(A33) 
(A34) 
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